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Abstract 

Using the Feigenbaum renormalization group (RG) transformation we 
work out exactly the dynamics and the sensitivity to initial conditions for 
unimodal maps of nonhnearity > 1 at both their pitchfork and tangent 
bifurcations. These functions have the form of g-exponentials as proposed 
in Tsalhs' generahzation of statistical mechanics. We determine the q- 
indices that characterize these universality classes and perform for the 
first time the calculation of the g-generalized Lyapunov coefficient A,. The 
pitchfork and the left-hand side of the tangent bifurcations display weak 
insensitivity to initial conditions, while the right-hand side of the tangent 
bifurcations presents a 'super-strong' (faster than exponential) sensitivity 
to initial conditions. We corroborate our analytical results with a priori 
numerical calculations. 

PACS numbers: 05.10.Cc, 05.45.Ac, 05.90.+m 

The nonextensive generalization [1] of the canonical statistical mechanics has 
raised interest in testing its applicability in several suggested circumstances in 
a variety of physical systems [1]. A class of problems where this issue has been 
much studied recently is the dynamical behavior of nonlinear iterated maps un- 
der critical conditions [2, 3, 4, 5, 6]. Such is the case of the bifurcation points 
associated to deterministic chaos in simple nonlinear dissipative maps, like those 
occurring in the logistic map and its generalization to nonlinearity C > 1- These 
types of critical states provide an exceptional opportunity to examine explic- 
itly the underlying mathematical structure and physical implications of their 
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Figure 1: Schematic form of at the a) pitchfork and b) tangent bifurcations. 

scale-invariant, power-law properties. As it is well-known the period-doubling 
and intermittency routes to chaos are based on the pitchfork and tangent bifur- 
cations, respectively [7]. The pitchfork bifurcations are the mechanism for the 
successive doubling of periods of stable orbits, and their accumulation point, an 
orbit of infinite period, is the so-called chaos threshold. The tangent bifurcation 
is a different mechanism linking periodic orbits with chaos in which intermit- 
tency is a precursor to periodic behavior [7]. Like in conventional equilibrium 
second order transitions, the bifurcation points have universal properties that 
can be obtained by means of renormalization group (RG) techniques, which 
for maps of the logistic type are based on functional composition and rescal- 
ing. As a result of this the static universal properties such as those derived 
from fixed-point maps and their perturbation equations are well understood [7] . 
However, the dynamical properties and their implications are only now being 
studied in a similar fashion by the same RG approach [8]. Here we present RG 
analytical results for the universal dynamics at bifurcation points of unimodal 
(unidimensional with a single maximum) maps of arbitrary nonlinearity C > 1- 
At each of the map bifurcation points the sensitivity to initial conditions 
obeys power law behavior for large iteration time t in contrast to the usual 
exponential growth or decay generic of noncritical states [2] . This of course im- 
plies the vanishing of the Lyapunov Ai coefficient that measures either periodic 
(Ai < 0) or chaotic behavior (Ai > 0). The generalized nonextensive theory 
offers an attractive alternative to describe the dynamics at such critical points 
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[2], the sensitivity would in these cases obey the q-exponential expression 
with a g-generaUzed Lyapunov coefficient Xq 



^t=expg{Xgt) = [l-{q-l)Xgtr— (geR), (1) 

that yields the customary exponential with Ai when g — > 1. The definition 
of Xq differs from other generalizations of a Lyapunov coefficient [9]. One of 
us has recently pointed out [8] a clear connection between the fixed-point RG 
solutions at bifurcation points of unimodal maps and nonextensive entropy ex- 
tremal properties. The aim of this letter is to work out analytical solutions for 
the dynamics and the sensitivity to initial conditions in the proximity of the 
first bifurcations of both types, and to corroborate these results numerically. In 
doing this we confirm previous estimations of the nonextensive parameter q and 
calculate for the first time the g-generalized Lyapunov coefficient A, in terms of 
the map parameters. 

Let us start by briefly recalling the background formalism. We consider the 
generalization to nonlinearity of order ^ > 1 of the logistic map 

f^ix) = 1 - (2) 

where — 1 < a; < 1 and < /x < 2. These maps display a "pitchfork bifurcation 
regime" for /x < Hoc where /ioo is the value of the period-doubling accumulation 
point (onset of chaos). On the other hand, as a consequence of the tangent 
bifurcations the "chaotic regime" for jj, > /Xoo becomes interrupted at certain 
values of by windows of periodic behavior [7] . The solution of Hu and Rudnick 
[7, 10] to the Feigenbaum RG recursion relation [11] for the tangent bifurcation 
was obtained as follows. For the transition to a periodic window of order n 
there are n points for which the original map / is tangent to the line of unit 
slope. Choosing one of these points, shifting the origin of coordinates to this 
point, and making an expansion of the n-th composition of the original map, 
one obtains 

&\x)=x + u\x\' + o{\xn, (3) 

where u is the leading expansion coefficient. The RG fixed-point map x' = f*{x) 
was found to be 

x' = x[l -{z- l)u sgn{x)\x\^-^]-^ . (4) 

This solution has a power-series expansion in x that coincides with Eq. (3) in 
the two lowest-order terms. In Ref. [8] it was observed that the previous scheme 

is also applicable to the pitchfork bifurcations rf/ ^\x)/dx = — 1 of order 

x=0 

n = 2*^, fc = 1, 2, provided that the sign of u is changed for a; > 0. As it will 

become clciar below, wc note that the power z appearing in Eqs. (3) and (4) is 
different from the power ( in the original map Eq. (2). 
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Now, taking into account that the fixed-point map Eq. (4) satisfies /* (/* (a;)) = 
a~^f*{ax) with a = 21/(^-1), we obtain by repeated functional composition the 
foUowing remarkable property 

r'"'\x) = ^rim^^x), m = l,2,... (5) 

This property implies that, for a total number of iterations t = mn, m = 1,2, 
with sufficiently small initial Xq, and for n and therefore t large, the fixed-point 
map can be written as 



■ {xo)=xo[l-{z-l)asgn{xo)\xorH]-^, (6) 



(m) 



where a = u/n. This is equivalent to a continuous time t approximation of Eq. 
(4) obtained by taking xt = x' , xq = x and at = u. We notice also that by 
introducing an additional index I = l,2,...,n — 1, the same expression holds 
for Xij^i = /^™"~'''-'(a;o), since .tq on the right-hand side of this equation can be 
replaced by f^''\xo)- In words, the time evolution of the original map in the 
neighborhood of a bifurcation point of order n, consists of a bundle of n orbits 
or trajectories of the form (6). But if one considers instead the time evolution of 
the n-composed map /^"^ this is given by a single orbit of the form (6), only now 
with t = m and a = u. It can also be verified that Eq. (6) satisfies the property 
dxt/dxQ = [xt/xoY , and this implies that the sensitivity to initial conditions 
= limAa:o^o(^a;t/Aa;o) has the form 

Uxo) = [1 - - l)a sgn{xo)\xo\'-H\-^^ . (7) 

Comparison of Eq. (7) with Eq. (1) yields the identifications [12] 

q = 2-- and Aq(a;o) = zo sgn(a;o)|a;o|''~^ (8) 
z 

Further, by using the known [13] form p{x) ^ \x\ for the invariant distri- 

bution of /'"^ in Eq. (3) , we have that the average Ag of Xq{xQ) over xq yields 
the piece-wise constant 

Xq = za sgn(a;o). (9) 

It is interesting to notice [8] that this average corresponds to the (/-extension 
of the customary expression for the Lyapunov coefficient Ai, obtained as the 
average of In \df{x)/dx\ over p{x). Here we have 



Ag = dXp{x) lUg 



df''\x) 



dx 



(10) 



where lng(2/) = {y^ ^ — 1)/(1 — q), and In^ (exp^(x')) = x. Likewise, for the 
n-composed map /^"^ the sensitivity to initial conditions is characterized by the 
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Figure 2: Weak insensitivity to initial conditions for the 1st pitchfork bifurcation 
when C = 1.75. The straight hne for the g- logarithm of vs t with g = 5/3 
and the slope m — —2.55 confirms the analytical result and the specific values 
predicted for q and \q for this transition. The inset shows the corresponding 
behavior for a single trajectory. 



g-generalized Lyapunov coefficient A, = nXq 



In Ref. [8] the alternative 



expression for the sensitivity to initial conditions was proposed 
6 = [exp,(A,t)]' = [1 - (9 - l)A,i]-^, 



(11) 



as suggested by a noteworthy similarity between the RG perturbation expression 
for xt and the corresponding nonextensive generalization of the Kolmogorov- 
Sinai entropy [8]. Using this form, the identifications are q = z and \q = a. 

Next, we proceed to show explicit results for the index q and the g-generalized 
Lyapunov coefficient Xq for the first pitchfork and tangent bifurcations. A pitch- 
fork bifurcation of order n = 2^ (see Fig. la) is characterized by the following 
conditions 



dy 



y=Vc 



Vc (fc = l,2,...). 



(12) 



— — 1 and 



dy 



V=Vc 



these determine the values of the critical parameter fic and of the 2 critical 
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positions yc- Shift of coordinates [x = y — Uc and x' = /^^ \y) — yc) and 
expansion yields 

/^^"^ (x) =x + u\x\^ + o(a;^), (13) 

since cP' f^'^''^ / dx'^\x=ci = always at these transitions, independently of the 
nonlincarity The coefficient u is then given by u = ±{l/6)d^f'-^'°^/dx^U^o, 
where the + sign applies for a; > and the — sign for a; < 0, and where 
d^fC^ ^/dx^\x=o < 0. As a first result we have that z = 3 and q = 5/3 for the 
sensitivity to initial conditions of all the pitchfork bifurcations of any nonlincar- 
ity Note that for any value of xq ^ the argument of the g-exponential (7) is 
negative, and this, combined with g' > 1, implies that the pitchfork bifurcations 
display weak insensitivity [2] to initial conditions at both sides of the critical po- 
sition. In the case of the first pitchfork bifurcation [k = 1) its location is given 
by /ic = (C + 1)''/(C''(C + 1)) and yc = C/(C + 1)) whereas the coefficient u is 
given by w = =f(C^ + 2^^ ~ 2^ — 1)/(6C^). So that the ^-dependent ^-generalized 
Lyapunov coefficient for /(^) is 



In Fig. 2 we present the numerical results for the first pitchfork bifurcation 
with ( = 1.75 that corroborate the analytical prediction Ag = —2.5466. We 
have verified Eq. (14) for several values of C,. 

The location and the expansion coefficient for the second pitchfork bifur- 
cation {k = 2) for general are harder to obtain analytically, but as an il- 
lustration we give results for the logistic map (( = 2). We obtain /ic = 5/4 
and the two critical positions yc = 2(1 ± •\/2)/5. Interestingly, Xq for the first 
critical position is = —3(250 -|- 125\/2)/4, while for the second position is 
Xg = -3(250- 125^/2)/4. 

The conditions that determine the critical parameter and the n critical po- 
sitions of a tangent bifurcation of order n (see Fig. lb) are 



f-(ra)/.. ^ _ ] "J 



r\yc) = yc and 



dy 



= 1, (15) 

y=Vc 



so that, by carrying out the same change in coordinates as before, one obtains 
the expansion 

/(")(s) =x + u|a;|2 + o(a;2), (16) 

where the coefficient u is now given by u = {l/2)d'^ f'^"'^ /dx^\x=o > 0. Hence we 
have 2 = 2 and q = 3/2 for the sensitivity to initial conditions of all the tangent 
bifurcations regardless of the value of (. The left-hand side of the tangent 
bifurcation points exhibits, in analogy to the previous case, a weak insensitivity 
to initial conditions. However at the right-hand side of the bifurcation the 
argument of the g-exponential becomes positive and this, together with q > 1, 
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Figure 3: a): Weak insensitivity to initial conditions at the left-liand side of the 1st 
tangent bifurcation when C, = 2. The straight line for the g-logarithm of vs t with 
q = 3/2 and the slope m = —31.15 confirms the analytical result and the specific 
values predicted for q and A, for this transition. The inset shows the corresponding 
behavior for a single trajectory, b): 'Super strong' sensitivity to initial conditions at 
the right-hand side of the 1st tangent bifurcation when C, =2. The straight line for the 
g-logarithm of vs t with q = 3/2 and the slope m = 31.22 confirms the analytical 
result and the specific values predicted for q and \q for this transition. The inset shows 
the logarithm of vs t where faster than exponential growth can be clearly observed. 
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results in a 'super strong' sensitivity to initial conditions, i.e. a sensitivity that 
is faster than exponential. Once again we use the logistic map C = 2 as an 
illustration, and wc choose for its first tangent bifurcation (n = 3) at /ic = 7/4 
one of the three critical positions, yc = 0.031405... , and determine its expansion 
coefficient to be m = 15.608... . Thus for /^^^ the g'-generalized Lyapunov 
coefficient is Xq — ±31.216... for .xq ^ 0. In Fig. 3a) we plot the numerical 
results for the left side of this bifurcation, while in Fig. 3b) we show the 'super 
strong' sensitivity to initial conditions characteristic of the right side of the 
tangent bifurcation. 

Naturally, only for very low order n of the bifurcation and/or special values 
of the nonlinearity C it is feasible to obtain algebraic solutions for the locations 
of the transitions and then for the leading expansion coefficient u. Nevertheless, 
it is always possible to implement numerical methods to determine them, and 
consequently the g-generalized Lyapunov coefficient A,. As shown in Figs. 2 
and 3 the a priori numerical calculations for the simplest examples selected here 
provide a striking confirmation of our RG predictions for both q and Ag. 

In summary, we have fully determined the dynamical behavior at the pitch- 
fork and tangent bifurcations of unimodal maps of arbitrary nonlinearity C > 1- 
This was accomplished via the consideration of the solution to the Feigenbaum 
RG recursion relation for these types of critical points. Our study has made use 
of the specific form of the ^-logistic map but the results have a universal validity 
as conveyed by the RG approach. The RG solutions are exact and have the an- 
alytical form of g'-exponentials, we have shown that they are the time (iteration 
number) counterpart of the static fixed-point map expression found by Hu and 
Rudnick for the tangent bifurcations and that is applicable also to the pitchfork 
bifurcations [8] . The g-exponential form of the time evolution implies an analyt- 
ical validation of the expression for the sensitivity to initial conditions suggested 
by the nonextensive statistical mechanics [2]. It also provides straightforward 
predictions for q and A^ in terms of the map critical point properties [8] . We 
found that the index q is independent of and takes one of two possible values 
according to whether the transition is of the pitchfork or the tangent type. The 
g-generalized Lyapunov coefficient A^ is simply identified with the leading ex- 
pansion coefficient u. As wc have shown thc^sc predictions are unquestionably 
corroborated by a priori numerical calculations. Significantly, both families of 
bifurcations display non-canonical dynamical behavior (either weak insensitiv- 
ity or 'super strong' sensitivity to initial conditions). The universal attribute of 
our results invites some reflection on the rationale for non-canonical dynamics. 
This is most apparently linked to the fact that the tangency shape of the map 
at these critical points either effectively confines or expels trajectories causing 
an abnormal, incomplete, sampling of phase space (here — 1 < x < 1). 
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